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I. INTRODUCTION 



High-energy scattering cross sections are usually expressed in a factorized form, in which 
the hard scattering kernel and the parton distribution functions (and the fragmentation 
functions if semi-inclusive or exclusive scattering is considered) are convoluted. The factor- 
ization property reflects the idea of separating long-distance and short- distance physics in 
high-energy scattering. The hard scattering occurs independent of the details about how 
energetic partons are produced or how final-state particles hadronize, and vice versa. The- 
oretically the factorization implies that each factorized part can be treated separately. The 
hard scattering kernel can be computed using perturbative QCD. Parton distribution func- 
tions describe the nonperturbative nature of partons to have a certain fraction of longitudinal 
momentum inside an incoming hadron, hence not computable in perturbation theory. But 
their evolution with respect to the factorization scale \x can be computed. The factorization 
property is the basis to study many high energy processes. 

In inclusive high-energy processes, only integrated parton distribution functions (i.e., 
parton distribution functions integrated over the transverse momentum) are needed to be 
convoluted with the hard scattering kernel. The integrated parton distribution functions can 
be well defined in a gauge-invariant way as the matrix elements of the relevant operators 
and their evolutions are described by the renormalization group equation. If we consider 
semi-inclusive or exclusive processes, that is, if one or more final-state particles are tagged, 
more detailed information on the scattering such as the transverse momentum dependence 
or spin correlations can be probed, but it is more complicated than inclusive processes. In 
this case, the scattering cross section is convoluted with one more nonperturbative frag- 
mentation function, which contains the information on the hadronization effect. However, 
since the transverse-momentum-dependent (TMD) parton distribution function (PDF) was 
proposed Q] , its definition and, accordingly, its radiative corrections are known to cause 
severe problems [Jij]. These problems are related to the failure of showing the cancellation of 
divergences between real gluon emission and virtual gluon exchange in certain calculational 
schemes once the TMD PDF is defined. The issues can be categorized into the definition 
of the TMD PDF, the choice of gauge (axial or covariant), and the regularization schemes 
though they are intertwined with each other. In this paper, the aspect of the transverse 
momentum dependence in semi-inclusive deep inelastic scattering (SIDIS) is considered and 
the spin correlation is not treated. TMD PDF has attracted a lot of attention since its effect 
can be detected in experiments. 

The problem can be described as follows [3]: The TMD (or unintegrated) quark distri- 
bution function f(x, kj_) at one loop is written as 



where x is the lightcone momentum fraction and k^ is the transverse momentum. The 
quantity P±r is the real gluon emission amplitude at order a s , schematically given by 



/« (x, kj.) = P 1R (x, kj.) - 5(1 - x)5(k ± ) / dx'd 2 k' ± P lR (x', k' ± ), 



(1) 



Pir{x, k±) ~ a s 



r l l 



+ regular as x — » 1 , 



(2) 



Ll -xk 2 , + A 2 
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where A 2 involves the infrared regulator such as small quark and gluon masses, and the 
second term in Eq. corresponds to the virtual gluon emission. It is argued that for any 
kj_, there is an endpoint singularity as x — > 1, and there is no cancellation between the 
real gluon emission and the virtual gluon contribution. If the calculation is performed in 
light-cone gauge, the divergence arises from the l/k+ singularity in the gluon propagator. 

There have been several approaches to solving this problem. First, Collins and Soper 
[l| proposed a non-gauge-invariant definition for the TMD PDF, and chose the axial gauge 
n ■ A = 0, but with a non-lightlike gauge fixing vector n M . As a result, a rapidity parameter 
( = (p ■ n) 2 /n 2 (n 2 ^ 0) was introduced and an additional evolution equation of the PDF 
on this parameter was developed requiring that physical quantities like scattering cross 
sections are independent of this scale. But the computation off the lightcone is complicated, 
and the factorization property and the relation to the integrated PDF are not clear in this 
approach (if. Collins and Hautmann 0] suggested to put Wilson lines on the lightcone and 
redefine the TMD PDF by adding additional soft Wilson lines to remove the divergence. 
They used method of regions to determine which kinematic regions yield divergences, and 
found subtraction terms to remove the divergence. From the explicit one-loop computation, 
they found additional soft Wilson lines slightly off the lightcone to be inserted in the gauge 
invariant definition of the TMD PDF. Cherednikov and Stefanis [fl redefined the TMD PDF 
by inserting a transverse gauge link which provides an additional soft counterterm. This 
additional term has the effect of cancelling the unwanted divergence related to the cusped 
contour. All the approaches try to define TMD or unintegrated PDF in a gauge-invariant 
way, and put some additional Wilson lines off the lightcone (taking the limit to the lightcone, 
if possible) in order to remove the problematic singularity mentioned above with the extra 
Wilson lines. 

Soft-collinear effective theory (SCET) 0, 0, 0] is the theoretical framework which is 
appropriate for describing energetic particles. SCET has been applied successfully in B 



decays and in other high-energy processes [10(] such as deep inelastic scattering [11], [l2|, [13, 
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151 ] . It can be also applied to SIDIS to describe the TMD PDF. The factorization property 
of SIDIS cross sections can be proved in a transparent way since the decoupling of collinear 
and soft degrees of freedom can be made manifest in the formulation of SCET. The TMD 
PDF can be defined in a gauge-invariant way putting all the Wilson lines on the lightcone 
without introducing extra Wilson lines off the lightcone. The computational technique 
developed in SCET offers an improved regularization method to show that the problem of 
the divergence is solved. The divergence appearing in Eqs. (CQ) and ([2]) is the mixture of 
ultraviolet and infrared divergences, which should not be present for physical quantities. In 
order to show the cancellation of the divergence using dimensional regularization, care must 
be taken to extract 5(1 — x) in Eq. (OQ). In this paper, the cancellation of the mixed divergence 
is explicitly presented. In Section [III kinematics relevant to SIDIS is briefly considered and 
the current operators in SCET for SIDIS are introduced. In Section [TTT| the TMD PDF and 
TMD soft Wilson lines are defined in terms of the matrix elements of the relevant operators. 
And the factorization of the hadronic tensor as a convolution of the fragmentation function 
and the TMD PDF is established. The TMD PDF itself is a convolution of the matrix 
elements for the collinear operator and the soft Wilson line. And the relation between the 
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TMD and integrated PDF is explained. In Section HVl the radiative corrections for the TMD 
collinear operator and the TMD soft Wilson lines are computed at one loop. In Sec. |VJ the 
renormalization group equations for the TMD collinear operator, the TMD soft Wilson line, 
hence for the TMD PDF are presented. In Section IVfl a conclusion is given. 



II. KINEMATICS 

In order to discuss TMD PDF, a reference frame should be selected first to define what 
transverse momentum is referred to. In SIDIS such as e + j) — > e + h + X , where h is 
an energetic tagged final-state hadron and X denotes all the remaining final states, the 
reference frame is chosen such that the spatial components of the momentum transfer g M 
from the leptonic system and the proton momentum P M lie in the z axis. In this frame, the 
momenta and P M can be written as 

(f = {n • q, q±, n-q) = (Q, T , -Q), 

P» = P£ = (n-P, P ± , n ■ P) ~ (A 2 /Q, T , Q), (3) 

where q 2 = —Q 2 , and A is a typical hadronic scale. The lightlike vectors n M and n M satisfy 
n 2 = n 2 = 0, and n-n = 2. We can set n M = (1, 0, 0, 1) and W = (1, 0, 0, —1), and the proton 
moves in the n (—z) direction with this choice. Note that an individual parton inside a proton 
can have nonzero transverse momentum of order A due to the fluctuation inside a proton, 
though the proton has no transverse momentum. We will consider the case in which the 
transverse momentum of the final-state particle is of order A. Of course, the outgoing hadron 
can have transverse momentum much larger than A. But in this case, we know from the 
momentum conservation that there should be other final-state particles with large transverse 
momenta such that the final transverse momenta add up to be of order A. Technically, as 
will be seen later, if the final hadron has a large transverse momentum, it is extracted as 
a label momentum. Then the large transverse momentum and the transverse momentum 
of order A are treated separately. Due to momentum conservation in each subspace of the 
transverse momenta, the final result can be derived in a similar way as in the case in which 
the final-state hadron is in the n direction with small transverse momentum of order A. For 
this reason, we consider the final-state hadron with n-collinear momentum P n and small 
transverse momentum from now on for simplicity, noting that the case with large transverse 
momentum can be treated in a straightforward way. 

At the hadronic level, we can define the invariants in terms of Q, P n and the outgoing 
hadron momentum P n as 

Q Q Pfi ' Pn ' Pn i a \ 

2P n ■ q n- P n P n -q Q 

Defining the momentum fractions £ and £' of the partons to hadrons as n ■ p n = £n ■ P n and 
£'n • p n = n-P n , the corresponding partonic invariants are given by 

_ x H Q _ z H _ n-p n 
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where the momenta of the incoming and outgoing partons are in the n M , directions 
respectively such that 



Pn = {n-pn,Pn±,n-p n ) ~ [-q-' A 



Q x. 

Pn = ( n • Pn, Pn±, n • Pn) ~ A, —J . (6) 

The scale Q sets the large scale in the system, and there is a hierarchy of scales in the 
momentum components such that 

p£~Q(A 2 ,A,l), p£~Q(l,A,A 2 ), (7) 

where A ~ A/Q, and SCET describes the dynamics of these particles. With this choice 
of the reference frame, the incoming partons and the outgoing parton can have transverse 
momenta of order A with respect to the axis on which the incoming proton and the photon 
lie. 

The electroproduction in SIDIS is mediated by the electromagnetic current operator of 
the form ifj'j^ip. The corresponding electromagnetic current operator in SCET is given by 

j^x) = c(g, / /) e i ^^-/ 2 -^"-/ 2 )x n n t 7±^Xn(^), 

jI(x) = C(Q,/i)e- J ^"-^ 2 -"^^ 2 )xn^ t 7±uV;x„(x), (8) 

where C(Q,fi) is the Wilson coefficient obtained by matching the current operator in the 
full theory onto SCET 
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Cm^) = l + ^(-ln 2 ^ + 31n^-8 + ^). (9) 
47r V n l n z 6 / 

A collinear fermion field in the n direction can be obtained from the full theory as 

a*) = X>-^fv^), (io) 

where = n ■ is the large label momentum. The field Xn, which is defined as 

Xn = W,l£n, is introduced to simplify notation, and W n is the collinear Wilson line given by 

W n = ^ eX P ~9=^ n ■ A n , (11) 



perm 



where n ■ V is the operator extracting the label momentum. For n collinear particles, the 
corresponding quantities are obtained by switching n and n. 

The collinear fields can be redefined to decouple from the soft fields as jsl 

£n(^) ► Y n ^ ni Af^ > Y n A^Yjl 

£,n(%) > Yn£,ni Aft ► Y^A^Y^ (12) 
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FIG. 1: The prescription for the soft Wilson lines in deep inelastic scattering, (a) In the current 
jjj,(x), £n from — oo to x, and £ s from x to oo. In j u (0), £,n from — oo to 0, and £ n from to — oo. 
See Eq. (I14j) . (b) The net soft Wilson line after the cancellation. 



where Y n and Y n are soft Wilson lines. The form and the analytic structure of the soft 
Wilson lines can be determined from the direction of the underlying collinear particles, and 
depend on whether the collinear particle is a particle or an antipartiele & And the bas.e 
building blocks for the soft Wilson lines are given as 



Y n (x) = Pexp 



dsn ■ AJns) 



Y n (x) = Pexp 



dsn ■ A s (ns) 



(13) 



and their hermitian conjugates, where P and P denote path and anti-path ordering respec- 
tively. In Eq. (jSJ), the soft Wilson lines are chosen appropriate for deep inelastic scattering, 
as explained in detail in Ref . [l6| . The prescription of the soft Wilson lines is shown in Fig. [T] 
(a), which is equal to Fig. [1] (b) since the overlapping Wilson lines are cancelled. 

There exist other equivalent prescriptions for the soft Wilson lines, but the procedure of 
calculating anomalous dimensions should also be prescribed. The above prescription is for 
the forward Compton scattering amplitude in deep inelastic scattering. (See Fig. [2j) This 
prescription is simple in the sense that the same Feynman rules can be used anywhere in 
Feynman diagrams, and it makes the calculation of the anomalous dimension straightfor- 
ward. Note that the choice of the directions for the underlying collinear particles such as 
forward scattering amplitude applies only to the prescription for the soft Wilson lines, not 
to the actual processes. In other prescriptions, the same result can be reached, but by a 




FIG. 2: Feynman diagram for forward Compton scattering amplitude in deep inelastic scattering. 
The scattering cross section is proportional to the discontinuity of the intermediate states by the 
optical theorem. 
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different procedure. For example, if soft Wilson lines are prescribed with a matrix element 
squared for a scattering cross section instead of the time-ordered product for the forward 
scattering amplitude, different Feynman rules should be applied across the cut to produce 
the desired result [lj], or the cusp angle in computing the anomalous dimension should be 
inverted in the complex conjugate part at the final step fl7| . 



III. TMD PARTON DISTRIBUTION FUNCTIONS AND SOFT WILSON LINES 

Consider the hadronic tensor, which is defined as 

WW = E / d 4 xe^(P\ 3 l{x)\Xh)(Xh\UV)\P), (14) 

where the sum is over all the final states X except the tagged hadron h. The SIDIS cross 
section is proportional to the hadronic tensor W^ u . It can be written in terms of the SCET 
current operators. The exponential factors extracting the label momenta all cancel due to 
momentum conservation requiring x = z = 1 (or £ = xh, £' = Zh), arid it becomes 



(y^Y^ (j ±vX n) d k (0)\Pn) (15) 



= C 2 (Q) d'xJ^i^lXnYhx^Xnix^X^iXhlxX^YnXnmP^ 
x 

= C\Q) f rf 4 x^(P,|(x,7±4)X^ y ab (I X 0- (x)|X/l)(X/l| (^ 



x 



where the color indices (a, b, c, d) and the Dirac indices (i, j, k) are explicitly shown in 
the last expression. The state \P n ) denotes the proton state in SCET, which consists of 
n-collinear partons. 

Since there are no n-collinear particles in the final state, the fragmentation function 
D(zn,n- r, rj_) is defined as the average of the following matrix element over color and Dirac 
indices such that 



^E< l(lx»(«))j^>W(^(0))jO> = ^ j ^e-^D(z H ,n.r,r ± ), (16) 



where D(zn,n ■ r,r±) depends only on n ■ r and r± and the label momentum is given 
by n ■ p n = n ■ P n /zH. The definition of the fragmentation function also holds when the 
final-state hadron h has large transverse momentum after summing over all the final-state 
particles, though the partonic variables x and z can have values less than 1, and appropriate 
adjustments should be made accordingly. Note that, as in inclusive deep inelastic scattering, 
when the contribution of the final-state particles is all summed over, there is no dependence 
on the transverse momentum and it becomes the jet function Jp(n ■ r) [12], which depends 
only on n-r. Therefore the transverse-momentum dependence in the fragmentation function, 
and subsequently in the PDF appears only when we specify at least one final-state particle 
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with nonzero transverse momentum. From now on, the subscript in zh (and Xh) is dropped 
for simplicity with the understanding that the momentum fractions refer to the hadronic 
variables. 

Using the fragmentation function, the hadronic tensor W^ u = —g^ v W can be written, by 
integrating over n ■ r, as 

W = j duC 2 (u) j dn 2 ^ ± j dn ■ rd 2 r ± e- t{n - rJl - x / 2+r ^D(z,n ■ r,r ± ) 

x (P n \x n (n-x,x ± )8(u-V + )^x*(0)\Pn)^ (17) 

where the soft Wilson lines are extracted from the n states, since they are decoupled from 
the collinear part, and the color projection is taken. In Eq. (fTTj) . the delta function 5{uj — V+) 
is inserted explicitly with V + = n ■ V + n ■ since the Wilson coefficient is actually an 
operator C(n-V,n-V^). By expanding the operators with respect to x±_ and n-x, we obtain 

X n (0,x ± )p(co-V + )xn(0) = J d 2 k ± e-^x^(uj-V + )5^(k ± + td ± )xn(0), 

tv (Tn . f <r , WiV-((W = / Anfie , „-i(vn-x/2+s ± -x ± ) 



itr Y^Y n (n ■ x, x ± )Y^Y n (0) = J drjd 2 s ± e 



x ltiYiY n 5(r] + n-td)5^(s ± + id ± )Y^Y n (0). (18) 

Note that the dependence of x n on n ■ s is dropped since the label momentum is already 
extracted, and the remaining momentum can be neglected. Then W is written as 

W = J duC 2 {u) J ^ J dn-rd 2 r ± J d 2 k ± J d V d 2 s ± 

= J duC 2 {uj) J dn-rd 2 r± J d 2 kj_ J dr]d 2 s±5(rj + n ■ r)5 {2) (r± + k± + Sj_) 

x D(z,n-r,r ± )(P n \O c (u,k ± )\P n )K(r],s L ), (19) 

where the TMD collinear operator O c (uj, k±) and the TMD soft Wilson line S(r], s±) are 
defined as 

O c (u,k ± ) = Xnp^-V + )5^(k ± + zd ± )xn, 
S(v,s ± ) = ^tTY^YJ( V + n-id)8^(s ± + id x )Y^Y fi , (20) 

and K{r),s ± ) = (0\S{r],sj_)\0). 

Note that the collinear particles and the soft Wilson line are not on the light cone x M = 
n M (n • x)/2 in Eq. (fTTj) . but are slightly off the light cone in the transverse direction x±. 
It is the expression with which all the approaches on TMD PDF agree, but the treatment 
of the TMD PDF and soft Wilson lines takes different paths from this point. Previously, 
these quantities were manipulated at this stage since it would presumably cause a serious 
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problem by putting the particles on the lightcone. Here all the particles are on the light 
cone, as in Eq. (I19p . and this is the starting point in this approach. It will be shown later 
that the radiative computations can be performed after putting all the collinear particles on 
the lightcone. That is, all the collinear particles are on the lightcone, so are the collinear 
and soft Wilson lines by expanding about the lightcone. As a result, Eq. ( 120|) is used for the 
TMD collinear operator, and the TMD soft Wilson line. 

The hadronic tensor can be expressed in terms of the TMD PDF. To define TMD PDF 
in SCET, let us recall the definition in the full theory, which is given as 

f(y,P±) = J dn ^^ e ^^ x ^^\P\^{n " x? x ^ w ^ n ■ x > x ^> ( 21 ) 

where ^[(n-a;, x±), 0] is the Wilson line connecting the two points to make a gauge-invariant 
operator, and y is the longitudinal momentum fraction of the parton n ■ p n = yn ■ P n . It 
is straightforward to write this definition in SCET by replacing the quark fields with the 
collinear fields and the Wilson lines with the corresponding collinear and soft Wilson lines. 
However, there is one more modification necessary, which is due to the way SCET describes 
disparate scales. Physically, the parton with the momentum fraction y can emit soft gluons 
before undergoing a hard collision, and the longitudinal momentum of the parton in a hard 
collision is given by yn ■ P n — k, where k is the longitudinal momentum of the soft gluons. 
In general, k is negligible since it is of order A, but it should be kept to guarantee the 
momentum conservation in each subspace of the momentum region. Note that the label 
momentum of order Q and the residual fluctuation of order A in SCET reside in different 
subspaces, and that is reflected, for example, in the treatment of the integral 

dze i{u-n- P )z e i{K-n)z = s ^ np I dze i(*-v)z = 8 u n . p (2^)8{K - n), (22) 



in which the integration over the label momenta of order Q (u and n ■ p) yields a Kronecker 
delta, but the integration over the residual momenta of order A (k and rj) yields a delta 
function. 

The TMD PDF in SCET is obtained from Eq. (ED) by switching to the fields in SCET, 
and by replacing yn ■ P n with yn ■ P n — k. As a result, the TMD PDF in SCET depends on 
an additional scale k, which describes soft gluon emission, and it is defined as 

f(y,pj_,K) = j "duo J ^0^ e K^yn-Pn+K)n- x /2-i Px . x± 

x (P R \x fi (x ± )6(u-V + )^Xn(0)\PK)±^ (23) 

where u = 2n ■ p n comes from the label momentum, and the soft Wilson lines are decoupled 
from the collinear sector. The dependence onn-xinxn is neglected because the remaining 
momentum of the collinear fermion in the n direction is much smaller than the soft mo- 
mentum in the soft Wilson lines. By expanding with respect to x±, f{y,p±,n) is written 
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as 

~} n~\ ,-vj ^7 2 



f{y,P±,K) = J du J dW 2 ' ^ g - L e »(^/2- y n-P a )T t - 3! /2 e i(^- a! /2+ Px . a!x ) 

x J d 2 k ± e- lk ^(P n \O c (tu,k ± )\P n ) J dr ] d 2 s ± e- l ^ x / 2 - iS ^K( V ,s ± ) 

= J du J d 2 k L J dr)d 2 s±5u,2yn-p n ${ K ~ v)^ (2 \p± + All + s±) 

x (P n \O e {u,k ± )\P fi )K(Ti,8 ± ), (24) 

where Eq. (I2"2l) is used. From the delta function 5(k —t]), it is clearly seen that k describes 
the soft gluon emission. To simplify further, let us define the matrix element of O c (u, k±) 
as 

(P n \Xnp(iU - V + )5^(k ± + id ± ) X n(0)\Pn) 

= j duj d 2 p ± s(u- ^p-) 5(2) (^ + P±Mu, P± ). (25) 

Plugging this expression into Eq. (I24p and integrating over the delta functions, the TMD 
PDF is a convolution of the matrix elements of the TMD collinear operator and soft Wilson 
lines, and is written as 

f(y,P±,K) = J d 2 p_ L d 2 s ± 5 ( - 2) (p ± - p ± + s±)g(y,p±)K(K,s±). (26) 

From the delta function, we can see that the transverse momentum of the parton p± is 
given by the original transverse momentum p±_ of the collinear parton before the soft gluon 
emission, subtracted by the transverse momentum of soft gluons s± {p±_ = p± — s±). From 
Eq. (125]) . the relation between g(y,p±) and O c (lu, k±) is given by 

g(y, = (Pn\O c (u = 2yn ■ P n , k±)\Pn). (27) 

In terms of the TMD PDF, the hadronic quantity W is written as 

W(y,z) = J duoC 2 {u) J dud(u—- — — ^ J dn-rd 2 rj_ J dr)5(n ■ r + r))D(z,n • r,r±) 



x 



/ d 2 s±d 2 p ± 5 ( - 2) (r± + s± - p±)g(u 7 p±)K(r], s±) 



where Eq. ff26l) and the fact that uo = In ■ p n = 2yn ■ P n are used in the last expression. W 
factorizes into the TMD fragmentation and the TMD PDF, and the TMD PDF is factorized 
in terms of a convolution of the collinear and the soft parts, as given in Eq. (1251) . Eqs. (|26|) 
and (I28p are the main result which proves the factorization of the hadronic tensor. In 
Eq. ( 1281) . the energy-momentum conservation is manifest from the delta functions. The 
small momentum components in the n direction in the hadron and in the soft gluons are 
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compensated by each other (n ■ r + rj = 0), and the transverse momentum p± from the 
incoming parton is equal to the sum of those for the outgoing hadron (r±) and the soft 
gluons (s±) with p± = + s±. 

It is useful to consider the relations between the TMD PDF and soft Wilson lines and the 
integrated PDF and soft Wilson lines, which appear in inclusive deep inelastic scattering. 
In SCET, the relation between the TMD PDF and the integrated PDF can be established 
in a transparent way. The integration of f{y,p±, k) over the transverse momentum can be 
performed easily when we consider the Fourier transform (or the impact parameter space 
representation) for g(y,p±) and K (rj, s±) 

g(y,b ± ) = [ d 2 p ± e^g(y,p ± ), g(y,p ± )= [ ^ e -^~g(y,b ± ), 



, (2tt) 2 ' 

K(k,c ± ) = J d 2 s ± e lc ^K(K,c ± ), K(k,s ± ) = J ^e- ic ^K(n,c ± ). (29) 
Then the impact parameter space representation of the TMD PDF is written as 

/(3/,&l,k) = J d 2 P± e ib ^f(y, P± ,K) = ~g(y,b ± )K(K,-b ± ). (30) 
By putting b± = in Eq. fl30l . we have 

f(y, K )= / d 2 p ± f(y,p ± ,K)= / d 2 p ± g(y,p ± ) / d 2 s ± K(K, s±) = g(y)K(n), (31) 



where the integrated collinear matrix element g(y) and the soft Wilson line K{rj) are given 
in SCET as [ljj 

M = j du5u,2yn-P n (Pn\Xf!^^ -V + )Xn\ P n) = J d 2 p±g(y,p ± ), 

K{rf) = ^(0\trY^Y n 6(r ] + n-id)Y^Y fl (0)\0) = J d 2 s ± K( V , s ± ). (32) 

It means that the integrated PDF is obtained by the product of the integrated collinear part 
and the integrated soft Wilson line to all orders in a s . The explicit relation at one loop 
between TMD and integrated PDF and soft Wilson lines will be presented in Sec. [V] 

Since the relation between the TMD PDF and the TMD collinear and soft Wilson line 
operators is established, we can focus on each operator and compute its radiative corrections. 
Note that the soft Wilson lines are present almost everywhere in high-energy processes 
though their analytic structure differs depending on the physical processes. Furthermore, 
as can be seen in Eq. (I15j) . part of the soft Wilson line (Y n ) comes from the collinear fields 
contributing to the fragmentation functions. Therefore it is hard to include the soft Wilson 
line residing only in the PDF. Instead they are sometimes regarded as universal objects, 
and are pulled out of the collinear part. And the collinear part only is referred to as 
contributing to the PDF. However, to make the comparison with the integrated PDF easy, 
the soft Wilson line is included in the PDF, which is consistent with the definition in the full 
theory. Therefore the distinction between the PDF and the collinear matrix element should 
be carefully understood from the context in the literature. 
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IV. RADIATIVE CORRECTIONS 



The relevant collinear and soft Wilson line operators in considering the TMD PDF are 
given as 

S(v, s±) = ^tr YlY n 8(rt + n ■ td)5 {2) (s ± + iV ± )Y^Y n . (33) 

From now on, the transverse momentum is written as a 3-vector. The radiative corrections 
for these operators and their renormalization group equations will be considered. 



A. One loop correction of the collinear operator 

Let us consider the radiative corrections for the collinear operator 

O c (u,k ± ) = Xn^-V + )6( D - 2 \k ± + tV ± ) Xn - (34) 

The Feynman rules for the operator O c to order g is shown in Fig. [3] by expanding the 
collinear Wilson lines W n . The Feynman diagrams for the radiative correction at one loop is 
shown in Fig.Hl The computation is performed with dimensional regularization in D = 4— 2e 
dimensions for the ultraviolet divergence and the nonzero external p 2 acts as an infrared cut- 
off. Note that the two-dimensional delta function in Eq. (1341) becomes (D — 2)-dimensional 
delta function in dimensional regularization. We can also use the two-dimensional delta func- 
tion. While the above definition satisfies the rotational invariance in the (D — 2) transverse 
plane, the two-dimensional delta function breaks this symmetry. Though the symmetry is 
broken using the two-dimensional delta functions, the result for the radiative corrections can 



(a) 




= fl -S(uj-n-(p + p'))^ D - 2 \k ± + p 1 _) 



= 9T„? [S(u - n ■ (p + p 1 ) + n ■ ^'^'{kx + P J 

2 n ■ q 

+8(uj - n ■ (p + p>) - n ■ 9 )5( fl - 2 )(kx + p, + qj] 
fi, a 

FIG. 3: Feynman rules for the collinear operator O c (uj,k±) to order g. In (b), T a are SU(3) color 
generators, and the momentum of the collinear gluon is incoming. 
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FIG. 4: Feynman diagrams of the radiative correction for the collinear operator O c at one loop. 



be made to be the same. Therefore the definition of the transverse delta function does not 
have to be unique, but it is preferable to keep the symmetry of the theory in dimensional 
regularization. 

Diagrams (a) and (b) in Fig. H] yield 



M a = M b = -Hg 2 C F % 



d D l n-(l+p) 



2 J {2ti) d l 2 {l + pfn ■ I 
-5{u -J -In- l)b {D ~ 2 \\^ + k ± + p ± 



(35) 



with {J = n-p + n-p'. The prescription for the poles is such that +i0 is added to each factor 
in the denominator. The first integral in Eq. ( 1351) is simple because the delta functions do 
not involve the loop momentum, and it is given as 



h 



d u l n-(l+p) r\c(D-2)n \ 



167T 2 V jJ? 



-p 



_eUV ^f-^-l)^-^ (2 - 2e) (kx). 

e UV ^ e IR ' 



(36) 



Note that there appear infrared poles of erR even though nonzero p 2 is introduced as the 
infrared cutoff. The second integral is more complicated and it is given by 

After integrating the delta functions, the remaining integral can be performed using the 
contour integral in the complex n-l plane. Depending on the position of the poles, the integral 
consists of two parts, each of which is proportional to 6{uj' — uj)6{uj) or 0{uj — uj')0{—uj). The 
first one corresponds to the contribution from a quark, while the second one is the antiquark 
contribution. Extracting only the contribution from a quark, the integral becomes 



i u 6{u' - uo)9{uj) 

167T 3 Uj' UJ 1 — UJ 



00 

k± + — p± 

io 1 



+ 



uj{uj' — uj) f—p 



(38) 



We can interpret the characteristics of the integrals 1\ and I2 by tracing where each term 
arises. In I\, the delta functions do not include the loop momentum. It means that 1\ 
comes from the combinations \n = or x n = V^n by extracting a collinear gluon and 

attaching it to the fermion on the same side with respect to the delta functions. On the 
other hand, I2 comes from the configuration in which a gluon is attached to the fermion 
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on the opposite side with respect to the delta functions. In other words, l\ corresponds 
to the virtual correction and I2 to the real gluon emission in the full theory. In inclusive 
deep inelastic scattering, the divergence of the form l/(euv e m) cancels between these two 



contributions 12j. However, comparing Eqs. fl36l) and fl38|) . the divergence does not seem to 



cancel. In the on-shell limit p 2 = 0, as uj' — > uj and — > — p_i_, I2 diverges as 

u'-u(k ± + P± )^ (39) 

but it does not cancel the divergence of the virtual gluon correction. Therefore it was claimed 
that if all the collinear particles were on the lightcone and the dimensional regularization 
was used, the divergence from the virtual and the real gluon corrections would not cancel 
unless an improved regularization technique was employed. 

Here an improved renormalization technique is suggested using dimensional regulariza- 
tion, in which the divergences between the virtual and the real gluon corrections cancel. The 
point is that care must be taken in extracting a delta function in dimensional regularization. 
In Eq. ( 138]) . the interesting limit is uj' — > uj and — > —p±. The naive limit in Eq. (|39|) 
does not extract the correct singular behavior in dimensional regularization. Note that there 
is another parameter — p 2 in Eq. (138j) . which is introduced as an infrared cutoff. This can 
be made arbitrarily small before any physical limits are taken. Therefore the basic idea is 
that we have to take the limit p 2 — > first, and then take other limits. The dimensional 
regularization changes the short- and long-distance behavior of the loop integral, and the 
transverse momentum and the n-component of the momentum are closely related to each 
other as a result. By putting p 2 = before any dependence on the dimension D = 4 — 2e is 
computed, it ruins this relation. In this respect, consider the limit 

lim t^—. (40) 
As it is, it diverges when kj_ — > 0. However, we can consider the function 

and take the limit a — > such that f(a, kj_) can be defined as a distribution function. As 
kj_ — > 0, /(a, kj_) ~ a~ 1+s and for kj_ 7^ 0, f(a, kj_) ~ a 5 . For < 8 < 1, f(a, kjj can be 
regarded as a representation of the delta function in the limit a — > 0, that is, 

\im-^ = A5^ D ~ 2 \k ± ). (42) 

The normalization constant A is independent of a, and is determined by the integral 

d °- 2k - I = J^L a ^ = (43) 
{27i) D - 2 k 2 ± + a (47T) 1 - (2vr) 2 - 2e ' 1 1 

If we choose 5 = e, A = 7r 1_IE r(e). Therefore as a — > 0, we can write 

^ 7r i-T(e)a-^( 2 - 2£ )(k ± ). (44) 

I ~\~ CL 
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This holds for < e < 1, and it can be analytically continued for other values of e. The pole 
is of the ultraviolet origin, as can be seen from Eq. (1431) . When the representation of the 



cannot 



delta function is used, it can be clearly seen that the momentum cutoff |kj_| < Q [18 
be used since the representation of a delta function is obtained only when all the range of the 
transverse momentum is included. Otherwise, there will be no cancellation of divergences, 
and the anomalous dimension will depend on the cutoff. Actually the integration should be 
performed for all in SCET, since the upper bound for the effective theory is pulled up 
to infinity, and the modification of the high-energy behavior is implemented in the Wilson 
coefficients. 

If Eq. ( |44"|) is applied to I 2 , it becomes 



i uj f—puj 
16n 2 uj' V [i 2 uj' 2 
-p 2 



"(u/ - u)- x - e 9(u/ - u)6{u))T{t)5 {2 - 2e) (k_L + -v± 

\ UJ 



(45) 



lQir 2 V /i 2 
We can write 



1 



-S(u — uj') + 



1 uj 6{J - uj) 
euv (a/ - lo) + 



9(uj)6^(k ± + - /P± 
x uj' 



(u/ - u)- 1 -^' - u)8(u) = (uj')- ^ - x) 



-l-e. 



'l-x)0{x) 



(46) 



where x = uo/uj' . When we integrate over x, the integral diverges at x ~ 1, which is of the 
infrared origin. Since the limits of the integral is bounded, there is no ultraviolet divergence 
and we obtain the formula 



{uj' - u)- l - e 9{uj' - uj)9{uj) 



1 , . 6(uj'-uj) 

— 5(u' -u) + -r 

e m (uj'-u}) + \ 



(47) 



with the -(—distribution function. In the above expression, the divergence of the form 
l/euv^iR is cancelled in the sum Ii + I 2 , and the amplitudes are given as 



M n = M h 



a s Cpi/i 2 



u) tf (w -aO^(kx) + ^f^^^(k ± + ^p ± ) . (48) 



4vr 2e uv 

In a similar way, M c can be calculated to give 



uj' (uj' — uo) + 



UJ 

up 



a s C F -2 



5 (2-2e)/ k 



UJ 



47r e uv 

Adding all the contributions, we have 

a s CF^ 1 



UJ 



UJ — UJ 



UJ' 



6{u' - lo)9{lo) 



(49) 



M a + M b + M c 



2n 2 euv 



+ ——, r V(U) - UJ 

(uj' -uj) + 



25{uo-uj')5 {2 - 2e \\i A _) 



(50) 



Now we have to subtract the zero- bin contribution 19j , which corresponds to the soft limit 
when the collinear loop momentum becomes soft, i.e., of order A. The collinear contribution 
from this soft region should be subtracted since collinear particles are defined to have large 
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label momentum. This region should be described by soft particles, not by collinear particles 
with small label momenta. The loop integrals in M a + M& + M c were performed naively 
including the contribution of the collinear particles with small label momenta. Therefore 
this contribution should be subtracted. In computing the naive collinear contribution, the 
loop and the external momenta scale as 

P = (n • M±, n ■ I) ~ (Q, A, A 2 /Q) , p" = (n • p, p ± , n ■ p) ~ (Q, A, A 2 /Q), (51) 

But in computing the zero-bin contribution, we subtract the contribution where the largest 
component n ■ I becomes A, retaining the same scaling for n ■ I and adjusting the scaling of l± 
such that n • In • I ~ l\. Therefore the power counting of the loop and the external momenta 
is given by 



F = {n-l,l ± ,n-l) ~ (A,Ay/A/Q, A 2 /Q), p» = (n-p,p ± ,n-p) ~ (Q,A,A 2 /Q). (52) 

One may wonder why the loop momentum of order I 2 ~ A 3 /Q is subtracted, hence unphys- 
ical. That is exactly the reason why it should be subtracted from the collinear contribution 
and replaced by the soft contribution in which all the components of the momentum are of 
order A. We can also choose the power counting as n ■ I ~ A, while others are fixed such 
that l\ 3> n • In • I, which is called the Glauber region. But it does not alter the result 
of the zero-bin contribution. The point of the zero-bin subtraction is to subtract the part 
in which energetic collinear particles become soft, and the criterion is the size of the label 
momentum. 

According to the power counting in Eq. (1521) . the zero-bin contribution from the naive 
collinear integrals M a and Mb are given as 

d D l 1 



x 



(2n) D l 2 (n-l+p 2 /n-p)n-l 
6(oo - u')6 {D - 2) (k ± + p _l) - 5(u -oj'-2n- l)5^ D - 2) (k ± + P± )l , (53) 



and the zero-bin contribution from M c vanishes. The first integral is given by 

16n 2 \fjLn- p) Veuv em^ euv 

The second integral J 2 can be computed by integrating over the delta function first, doing 
the contour integral in the complex n ■ I plane, and then finally integrating over l± using 
dimensional regularization. J 2 is given as 

% 1 6{uj'-uj) f(u' - lu)(-p 2 )\-^v 



nu-uA-vny™^ (k± + pj 



167T 2 euv W — UJ \ 71 • P/jL 

i ( -p 2 \-euv 1 r / 1 1 \., t\ 9{lo' - u)-\ -co 2f w, 
lb7r z \n-pfi/ euv L ^ euv em^ {uj' — uj) + \ 
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where the following formula is used: 
6(u' - w)(w' - u)- 1 -* 



euv 



uo) 



9{uo' 



00 



00 



(56) 



The difference between this equation and Eq. (H7|) is that there is no lower bound. Therefore 
Eq. (|56|) has both infrared and ultraviolet singularity. In the sum J\ + J 2 , the S(u' — uo) term 
cancels and the zero-bin contribution is given by 

a s CF^ 2 



2tx 2 euv \oo' 



(a/ - u) §(2 _ 2e) 



00 



(57) 



Subtracting the zero-bin contribution, the total collinear contribution with the wavefunc- 
tion renormalization is given by 



M a + M h + M c - (M° + M°) + w.f. 

asCpt/i 2 



2-n 2 euv 
1 + (u/uf) 2 



10J' 



00 



J -u)e{uj))5 (2 - 2e \k ± + p ± ) 



eu - w)^(cu)(5 (2 - 2e) (k ± + - P± 

to' 



(58) 



where 0(k>) in the second term is inserted for physical reasons. For comparison, the corre- 
sponding radiative correction for the collinear operator giving the integrated parton distri- 
bution function is given as [12[ 



agCpTfi 2 



2n 2 euv 



-6(u' 



uo) 



(uo' 



00 



U - uo)6(uo) + 1 + ^I^ QU - uj)6(uj) 



(59) 



This is the result from Eq. ( !58j) without the delta functions of the transverse momentum, or 
by integrating over the transverse momentum. However, note that the radiative correction 
for the collinear TMD operator is not just proportional to <5 ( - 2_2<E - ) (kj_ + p_i_), but it also 
depends on 5*- 2 ~ 2e ' (kj_ + (a//a>)pj_). This gives a nontrivial difference in the renormalization 
group behavior between the TMD and the integrated PDF. 



B. One loop correction to the soft Wilson line 

We now compute the radiative corrections for the TMD soft Wilson line operator 

S(v, s±) = ^tr Y^Y n 8( V + n ■ id)6® (s ± + iV x )Y^Y n . (60) 

The Feynman rule for the operator S with two gluons is given in Fig. El and the Feynman 
diagram for the radiative correction at one loop is shown in Fig. [61 The radiative correction 
is written as 

I( V ,s ± ) = -2ig 2 C F / -— ^ ; U K ±J 61 

J (271)^ U 2 (n ■ I - Ai)(n • I - A 2 ) 

6(rj + n ■ l)5^ D -V (s ± + 1 5{-q - n • Z)^" 2 ) (s ± - 1 ± ) 1 



(I 2 +i0)(n- 1 - Ai +i0) -n-Z + A 2 (P + i0)(-n • Z - Ai + iO) n • Z + A 2 J 
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v,b 



/i, a 



9 2N[{(n-q + it){n-q + it) + (n ■ q - it)(n ■ q + it) > [S± > 



q^ ->g 



nf'n" 



ni'rr" 



\(n ■ q + n ■ q — it) (— n ■ q + it)(n ■ q + it) 
+(a <-» 6, fi <-> j/, g <-> — g) 



FIG. 5: Feynman rules for the soft Wilson line operator S(r/, s±) with two soft gluons. 



where Ai and A2 are the infrared cutoffs. And the necessary iO prescriptions for the contour 
integral in the complex n ■ I plane are shown. 
The first integral is given by 

d D l 5(ri)5 {D -Vi 



Si 



1. / >■ 1 ' \ 1 1 n ' (Sj ) 

-^9 2 C F J (2n) D l 2 (n-l-\ 1 )(n-l-\ 2 ) 
a s C F 2 /AiA 2 ~ 



2tt efj y \ fi 



S{f])6 { 



2-2e), 



(62) 



In the second and the third integrals in Eq. (161]) . we perform the contour integration in 
the complex n ■ I place after integrating over the delta functions. The second and the third 
integrals are given by 

g 2 C F 0(17) 1 



^ 2 ^ 3 87T 3 77 + A 2 + 77A1 ' 
We apply the idea of the representation of the delta function again to have 



1 



lim 

Ai^O Si + 77A1 



7T 



1-e 



and the integral becomes 

a s C F 6{ri) (??Ai)- £ (2 _ 2e) 
27T euv ^ + A 2 



2vr 



A X A 2 \-«J_ 1 fl(ry) 



e uv e uv ^+ J 



5(2-2,), 



(63) 



(64) 



(65) 



Here the plus distribution function is also used. But note that all the poles in 1/e are of 
ultraviolet origin since the infrared cutoffs Ai and A 2 are used. Adding all these, the radiative 
correction for the soft Wilson line at one loop is given by 



J(r / ,s ± ) = ^^^ 2 --)(s ± ). 

7T e uv 7]+ 



(66) 




FIG. 6: Feynman diagram of one-loop correction for the TMD soft Wilson line. 
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The only difference, compared with the integrated soft Wilson line 16j, is the presence of 
8^ 2 ~ 2e '(s±). Therefore the anomalous dimension is the same as the integrated soft Wilson 
line except the delta function S^ 2 ~ 2e \s±), that is, the TMD soft Wilson line does not scale 
with respect to the transverse momentum. 



V. RENORMALIZATION GROUP BEHAVIOR OF THE TMD OPERATORS 



With the radiative corrections, the relations between the bare operators and the renor- 
malized operators can be written, in general, as 

0f(w,k L ;p_ L ) = J duj'd D - 2 k' ± Z c (uj,uj',k ± ,k' ± ;p ± )0 R (uj',k' ± ;p ± ), 

S b (v,sa_) = fdr/<p- 2 J ± Zs(ri,T/ t a 1 ,J ± )S R (T/M, (67) 

where the counterterm kernels Z c (uj, u', k±, k'j_; pj_) and Z s (r], r/, Sj_, s'j_) can be read off 
from Eqs. (1581) and (166|) . and are given as 



Z c (uj,uj',k ± ,k' ±] p ± ) = 8{u' -oo)8^\k x -k! ± 



+ 
+ 



2n 2 euv 
1 + (uj/uj') 2 
(u/ - uj) + 



(uj' — UJ)^ J 



k\ 



eas - u)9(u)5 {2 - 2e \k ± - k' - - — - P± ) 



Z 3 {ri,r/,aj_M = 8^(s ± - s' ± ) \s{ V - rf 



a s C F 9{r, - rf) 



(68) 



Here in Zc{uj, u', kj_, k'j_; pj_), the external transverse momentum pj_ of the quark is intro- 
duced for convenience, but it should be understood that it appears in the matrix element of 
the operator O c . The renormalization group equations for the renormalized operators are 
written as 

d 



VLj-0?(u,,1t ± ;p±) = ~ J ^W ±7c (c,u;',k ± ,k' ± ;p ± )Of(^,k' ± ), 
^S R ( V ,s ± ) = - J dri'd 2 s' ±ls (ri,ri', S± ,s' ± )S R ( V ',s' ± ), 

where the anomalous dimension kernels are given by 

a s C F 



(69) 



7c(u;,a; / ,kj.,k / _ L ;p±) 



7T 



-6(u-u/) K — J -^- 5 (2) k ± - k' ± 

2 (uj' — uj) + / 



+ \ + , {uj/uJ ' )2 6(cu' -uj)6(uj)5^(k ± -k' ± 



(u 1 - u) + 



UJ — LU 



ls{v,v', s ±, s ± 



71 



s± — s' 



(70) 



Note that the second term in 7^7 in Eq. (1701) . which corresponds to the zero-bin subtraction 
has the same form in 75 with the opposite sign. It confirms the nature of the zero-bin 
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subtraction, which is computed in the soft limit. This pattern also appears in integrated 
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It is useful to express the renormalization group equation for O c in terms of the dimen- 
sionless variables y and z defined as to = 2yn ■ P n , uj' = 2zn ■ P n , where < y, z < 1. It is 
given as 

^^Cc(y, k ± ; p ± ) = - J dz J d 2 k' ± 7c(y, z, k±, k' ± ; p±)O c (z, k' ± ; p ± ), (71) 



where 



lc{y, z, k ± , k' ± ; p ± ) 



a s C F 



+ y+wf svj^ -k' 



(z-v)- 
z-y 



P± 



(72) 



In terms of x = y/z, Eq. ( [71]) becomes 



djj, 7T 



+ 



1 + x 2 

:i-x) 



(73) 



Note that, from Eq. ff27j) . the matrix element of the collinear operator O c is given by 

g{y,-k ± ) = (P n \O c (cu = 2yn-P n ,k ± )\P n ). (74) 
Then Eq. ( I7ip can be expressed in impact parameter representation, or its Fourier transform 



as 



l*4-9(v,b±) = 1*4- I d 2 k ± e*>^g(y,-k ± ) 



d/j," 



a s C F 



(75) 



7T 



dz 



r3 



S(y-z) 



-2 + (1 + (y/zfje^-y/^^l _ 



0-y)- 



where the sign of the exponent becomes opposite compared to Eq. (|2"§|) since the 4-vector 
notation is changed to the 3- vector notation here. If we put = in Eq. (175|) . both sides 
become the integrated quantities over the transverse momentum k^, and the corresponding 
renormalization group equation is given by 



^IjlJ d2k±9 ^ y,k± " > = ~^ j ~( p gg( x ) ~ (i^) ) / d2 ^9{y/x,k } 



(76) 



where P qq {x) is the quark splitting function given by 



Pqq{ X ) — Cf 



3 f/i . 1+x 2 



(77) 
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An extra term in Eq. (ITS"]) other than P qq (x) comes from the zero-bin subtraction. Eq. (!76|) 
is exactly the same renormalization group equation for the integrated collinear operator in 
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The renormalization group equation for the TMD soft Wilson line after integrating over 



s'j_ in Eq. ( 169]) is written as 



d 



2a, C f 



7T 



drj 



iv - v')+ 



(78) 



Since the anomalous dimension for the TMD soft Wilson line is proportional to 5(s± — s 
the transverse-momentum-dependent part is trivial and the TMD soft Wilson line does not 
scale with respect to the transverse momentum. As a result, the renormalization group 
equation for the TMD soft Wilson line is the same as the integrated soft Wilson line 16 . 



Besides the renormalization group equation, the radiative corrections themselves, when 
integrated over the transverse momentum, give the same result for the radiative correction 
of the integrated collinear operator. If the counterterm kernels in Eq. ( 168]) is integrated 
over the transverse momentum, the results correspond to the counterterm kernels for the 
integrated collinear and soft Wilson line operators. Therefore the relation between the 
TMD (unintegrated) and integrated PDF can be clearly seen in this formulation. That is, 
the integrated PDF and its renormalization group behavior are obtained by integrating the 
TMD PDF and its renormalization group equation over the transverse momentum. We can 
see this relation explicitly. From Eq. ( 1261) . the renormalization group equation for f(y, p±, k) 
is written as 

/i^/(y, P_l, k) = J d 2 s ± d 2 k ± 5 {2) (p 1 _ + k ± + s±) 

x (v-^g(y, - k ±)) s -l) + 9(y, - k ±) (ji-jj-Kfa, s± < 
= - J d 2 s L d 2 k L 5 {2 \p ± + k ± + s±) J d 2 s' ± d 2 k' ± dzdn' 



x 7/(1/, z, k ± , k' ± ; k, k', s±, s' ± ; p±)g(z, -k' ± )K(K f , s 



(79) 



where 



7/(2/, z, k±, kj_; k, k', s±, s' ± ; p ± ) = 7 C (y, z, k ± , k' ± ; p±)5(k - /t')5 (2) (s± - s'J 

+ ls ( K , K ',s ± ,s' ± )5(y-z)5^(k ± -k' ± ). (80) 

This is a complicated integro-differential equation. But it is more convenient to derive the 
renormalization group equation in impact parameter space. From Eq. fl30l) . the renormal- 
ization group equation for f(y, b±, k) is given by 



^-^f(y,h ± ,K) 



^9(y, b±)) A-(k, -b ± ) + g(y, b ± ) (ji—K{k, -b ± ) ) . (81) 
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Each g(y,b±) and K(k, — bj_) satisfies the renormalization group equation 



dfx 
aJOx ' 



(82) 



H-^-K(k, -bj 
an 



71 

2a s C F 



dz 



-S(y-z) + 



-2 + (1 + (!//z) 2 )e 



2V<(1-1/A)bi-Pi 



7T 



2 ' ' (z-y). 

K — K' N 



dn'- —K(n', -b±). 

(k - «')+ 



Therefore Eq. ( !8~T1) is written as 
d[i 



dzdn ,; yf(y, z, k, b ± , p_l)/(z, bj_, «') 



(83) 



(84) 



where 



7 / (y,z,K,K / ;bj.,p J .) 



+ 



7T 

2^(k - «' 



[K — K 



(^(y-z 



-2 + (1 + (y/*)*)e 



2U«(l-!/A)bi-pi 



v«-2/)- 



)(5(/^-/^ , ) 
(85) 



This is the main result for the renormalization group behavior of the TMD collinear operator. 
For b^ = in Eq. flS5J), we obtain the renormalization group equation for the integrated 
PDF. 



VI. CONCLUSION 

It is shown that the hadronic tensor, hence the scattering cross section in SIDIS is fac- 
torized into a hard part, the TMD PDF, and the fragmentation function. The TMD PDF 
is further factorized into the TMD collinear part and the TMD soft Wilson lines. The 
TMD PDF is defined in terms of the matrix element of the gauge invariant TMD collinear 
operator in which all the collinear particles are put on the lightcone with no extra Wilson 
lines off the lightcone. The TMD soft Wilson is defined in a similar way. The radiative 
correction for the TMD PDF is computed with the zero-bin subtraction, and is shown that 
the divergences of the form l/(euv e m) are cancelled. This cancellation is obtained by using 
the representation of a delta function before the parameter e in dimensional regularization 
is set to zero. The resultant anomalous dimension kernel shows a nontrivial dependence 
on the momentum fraction x and the transverse momentum of the incoming parton. 
The radiative correction for the TMD soft Wilson line is also performed and the anomalous 
dimension kernel is the same as the integrated soft Wilson line except the delta function for 
the transverse momentum. That is, the TMD soft Wilson line does not scale with respect to 
the transverse momentum and it shows the same scaling behavior as that of the integrated 
soft Wilson line. 

There are several issues in obtaining the TMD PDF and its renormalization group be- 
havior. In contrast to previous approaches, all the collinear particles, hence all the collinear 
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and soft Wilson lines are on the lightcone. The cancellation of the divergence with the form 
1 / (euv e m) m the radiative correction of the TMD quantities is explicitly shown by carefully 
taking care of the limiting behavior of the calculation in dimensional regularization. Since 
the transverse momentum and the n-component of the collinear momentum are closely re- 
lated to each other in dimensional regularization, if we put the parameter e in dimensional 
regularization to zero in the first place, the calculation becomes inconsistent and the cancel- 
lation does not occur, as claimed before. The technique used here is to introduce external 
p 2 as an infrared cutoff and treat the ultraviolet divergence in dimensional regularization. 
Though there is an infrared cutoff, not all the infrared divergences are expressed in terms 
of the cutoff, but there are some divergences which appear as poles of epR. It is important 
to disentangle all the types of these divergences along with the ultraviolet divergence. After 
adding all the contributions, the infrared poles of eiR are cancelled. Especially, the problem- 
atic divergence is of the form l/(euv e m) which should not be present in order for the theory 
to make sense. In this approach there is no need to introduce additional soft Wilson lines 
which offers a counterterm to cancel the divergence from the virtual correction. Instead the 
use of the representation of a delta function and the plus function prescription are employed 
to show the cancellation of the unwanted divergence. 

In evaluating the radiative corrections, the external nonzero p 2 is introduced as an infrared 
cutoff and the limit p 2 — > is taken first before considering other limits. It is claimed that 
putting p 2 = before the calculation in dimensional regularization would yield an incorrect 
result. Then a question can be raised about whether the extraction of divergence is possible 
with the on-shell scheme (p 2 = 0) and treat both the infrared and the ultraviolet divergences 
using dimensional regularization. This is in principle possible, but it is extremely challenging 
to trace the origin of the divergences. In this case, however, the calculation with nonzero p 2 
can be a guideline to trace the divergences. 

The zero-bin subtraction is also important in avoiding double counting. In computing the 
radiative corrections of the collinear operator, the integrals were computed naively including 
the unphysical region with small label momentum. This region should be subtracted because 
that region is covered by soft particles. Otherwise, double counting occurs. The true soft 
contribution is decoupled from the collinear part, and it appears in the soft Wilson line. 

The renormalization group equation for the TMD collinear operator is a complicated 
integro-differential equation in which the momentum fraction x appears in the delta func- 
tion for the transverse momentum. In impact parameter space, the renormalization group 
equation involves an oscillating term which makes it difficult to solve the equation. However 
the renormalization group equations involve only physical variables. In contrast to previous 
approaches, this equation does not involve the rapidity parameter £ , or the mass param- 
eters m and A {ij, therefore it is difficult to compare the renormalization group equation or 
the radiative corrections directly. 

The comparison between the TMD (or unintegrated) PDF or soft Wilson lines and the 
integrated PDF or soft Wilson lines becomes manifest in this approach. In this paper, it 
is shown by explicit calculation at one loop that the radiative corrections and the renor- 
malization group equations for the integrated PDF and soft Wilson line are obtained by 
integrating those for the unintegrated (or TMD) PDF and soft Wilson lines over the trans- 
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verse momentum. This is different from previous approach in which the comparison was not 
transparent because the order of integrating over the transverse momentum and the removal 
of the divergence makes the comparison difficult ji)]. However the calculational procedure in 
this paper makes the comparison clear. 

Here the TMD PDF is considered only for quark distribution functions. But it is straight- 
forward to extend the idea, and calculate the radiative corrections for gluon distribution 
functions and their mixing at higher orders in a s . The analysis on the TMD collinear, soft 
operators and their renormalization group behavior can be extended to other high-energy 
processes such as jet production or Drell-Yan processes (for full QCD approach, see Refs. [l| 
and [20']). There will appear appropriate collinear and soft operators, but the form and 
structure of these operators might be different. It would be interesting to consider other 
high-energy processes and to see if there are important processes in which the transverse- 
momentum-dependent effects are relevant. 
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